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Abstract: We study the alignments of the low multipoles of CMB anisotropies with 
specific directions in the sky (i.e. the dipole, the north Ecliptic pole, the north Galactic 
pole and the north Super Galactic pole). Performing 10^ random extractions we have found 
that: 1) separately quadrupole and octupole are mildly orthogonal to the dipole but when 
they are considered together, in analogy to [7], we find an unlikely orthogonality at the 
level of 0.8% C.L.; 2) the multipole vectors associated to £ = 4 are unlikely aligned with the 
dipole at 99.1% C.L.; 3) the multipole vectors associated to £ = 5 are mildly orthogonal to 
the dipole but when we consider only maps that show exactly the same correlation among 
the multipoles as in the observed WMAP Syr ILC, these multipole vectors are unlikely 
orthogonal to the dipole at 99.7% C.L.. 
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1. Introduction 

The anisotropy pattern of the cosmic microwave background (CMB), obtained by Wilkinson 
Microwave Anisotropy Probe (WMAP), probes cosmological models with unprecedented 
precision (see [1] and references therein). WMAP data are largely consistent with the 
concordance A cold dark matter (ACDM) model, but there are some interesting deviations 
from it, in particular on the largest angular scales. Among these deviations, we focus on the 
alignments of low multipoles. Unlikely (for a statistically isotropic random field) alignment 
of the quadrupole and the octupole is described in references [2, 3, 4, 5, 6]. Moreover, both 
quadrupole and octupole align with the CMB dipole [5, 7]. Other unlikely alignments are 
described in [8]. 

We study directionality anomalies considering estimators defined for each multipole 
i. In the present paper we limit our analysis to the range A£ = 2 — 7. We take into 
account two different Monte-Carlo (henceforth MC) simulations. In the first one, we study 
directional anomalies considering the coefficients of the spherical harmonics expansions 
randomly extracted from a Gaussian distribution. In the second one we consider maps in 
which the relative alignments of the low multipoles are fixed and given by the WMAP 5yr 
ILC map (henceforth called ILC map for sake of brevity) . Both MCs are performed taking 
into account WMAPs anisotropic pixel noise (see below for details). Through the first MC, 
we obtain that quadrupole and octupole are mildly unlikely orthogonal to the dipole (this 
anomaly become stronger when they are combined in analogy to what has been performed 
in [7]). Moreover we find that £ = 4 is unlikely aligned with the dipole and ^ = 5 is very 
mildly unlikely orthogonal to the dipole. Through the second MC, we find that quadrupole 
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and octupolc are (mildly) unlikely orthogonal to the dipole and combining them together 
(still in analogy to [7]) we do not find any increase of this anomaly. In this case £ = 4 is 
unlikely aligned with the dipole and £ = 5 is unlikely orthogonal to the dipole. We do not 
obtain unlikely alignments with the dipole for ^ = 6 and 7 for both MCs. 

The paper is organized as follows: in Section 2 we present the Multipole Vectors 
expansion and the tools that have been used to perform the statistical analysis; such 
analysis is given in Section 3 where we present also our main results; conclusions are 
drawn in Section 4. Some technicalities about the performed simulations can be found in 
Appendix A. 

2. Tools and estimators 

The alignment of multipoles can be defined using a new representation of CMB anisotropy 
maps where the a^m (coefficients of the expansion over the basis of spherical harmonics) 
are replaced by vectors [2, 3]. In particular, each multipole order £ is represented by £ unit 
vectors Ui and one amplitude A 

aem^ A^^\ui, ■ ■ ■ ,ue. (2.1) 

Note that the number of independent objects is the same in the l.h.s and r.h.s. of equation 
(2.1): 2£ + l for equals S£ (numbers of components of the vectors) +1 (given by ^W) 
—£ (because there are £ constraints due to the normalization conditions of the vectors). 
For a more detailed explanation of equation (2.1) and of the properties of that association 
see for example references [2, 3, 8, 21]. 

Unfortunately, an explicit analytical expression for the association given in equation 
(2.1), is possible only for i = 1 and for £ ^ 1 numerical method are needed^. The Copi et 
al.'s algorithm (which use is acknowledged here) for constructing multipole vectors from a 
standard spherical harmonic decomposition is described in [2] and the implementation of 
it is publicly available^. Other methods exist but, as far as we know, their implementation 
is not publicly available on a standard platform (see for example [3, 9] where the problem 
of finding £ vectors is translated into the problem of finding the zeros of a polynomial of 
degree 2£). 

In order to investigate whether a map presents unlikely directions with respect to 
random CMB sky extractions^, we define the following estimators (scalars quantities) for 
each multipole order £: 

e 

Sei = Y,\^i-d\/£, (2.2) 

i=l 

where Ui represent the multipole vectors associated to a multipole order ^ of a map, d is a 
fixed direction in the sky and where the absolute value is taken into account since multipole 

^Indeed, ior £ = 2 it is possible to obtain the multipole vectors computing the eigenvectors of a symmetric 
and traceless tensor representing the quadrupole, see [11, 12]. 

^http:/ /www. phys.cwru.edu/projects/mpvectors/ 

^See Section 3.1 for a detailed description of the performed extractions. 
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vectors are vectors without "head" (this is the so called reflection symmetry) . The division 
by i is considered in order to define the estimators (2.2) between and 1. 

For sake of comparison, in this paragraph we present the estimator used in [5, 7] 
(where it is shown that quadrupole and octupole are unlikely aligned with the dipole). 
The definition of that estimator is the following 

where q and Oj are the area vectors built from the quadrupole and octupole multipole 
vectors respectively. More explicitly: 

^= 921 X 922 , (2.4) 

and 

o\ = 632 X 033 , (2.5) 

02 = 033 X 031 , (2.6) 

03 = 031 X 032 , (2.7) 

where q2i with i = 1,2 are the two multipole vectors associated to the quadrupole and 63^ 
with j = 1,2,3 are the three multipole vectors associated to the octupole. In eq. (2.3) the 
division by 4 has been taken into account to define the estimator 5321 in the range [0, 1] 
and the absolute values are considered to fulfill the reflection symmetry. 



3. The Analysis 

3.1 Description of the Random Extractions 

As already stated in Section 1 we perform two kinds of MCs in order to study how likely 
the ILC map ^ is aligned to some specific direction. 

For the first MC we extract 10^ random maps from a Gaussian distribution from an 

angular power spectrum corresponding to the best fit of WMAP Syr (see footnote 4). For 
each random extraction, the estimators S^i (with I = 2,3, 4, 5, 6, 7) and S'321 are computed 
for some direction of interest such as the Dipole (DIP), the north Ecliptic pole (NEP), the 
north Galactic pole (NOP) and the north Super Galactic pole (NSGP). 

For the second MC, as done in [5], we freeze the relative alignments of the low multi- 
poles such as those that are observed in the ILC map and perform 10^ random rotations of 
that map. For each random rotation the estimators Sn (with i = 2,3, 4, 5, 6, 7) and S'321 
are computed for the aforementioned directions (DIP, NEP, NGP and NSGP namely). 

Both MCs are performed taking into account WMAP's anisotropic pixel noise^. The 
resolution that has been considered for the analysis is given by the HEALPIx^ [10] param- 

*The 5 years WMAP ILC map, as well as other CMB data products, is publicly available at the Lambda 
web site: http://lambda.gsfc.nasa.gov/. 

^We consider the WMAP V band noise and hits files available at Lambda web site. See footnote 4. 
®http://healpix.jpl. nasa.gov 
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ctcr Nside = 16 that corresponds to 3072 pixels^. See appendix A. 2 for details about the 
degradation issue. 

3.2 Results 

In this Section we provide the results for both MCs. 

Results for the first MC are summarized in Table 1 where the probabilities to obtain 
a smaller value for the analyzed directions are listed for the ILC map. 

Table 1: First MC: Probabilities (in percentage) to obtain a smaller value for the considered 

directions and estimators 

Estimator DIP NEP NGP NSGP 



S'21 


6.7 


42.3 


9.6 


41.5 




3.6 


86.0 


7.2 


56.2 


anew 
"-"321 


0.8 


73.6 


2.6 


49.4 


S41 


99.1 


2.6 


48.2 


74.5 


S5I 


9.9 


30.9 


56.5 


87.7 




76.2 


58.8 


37.5 


95.0 


Sri 


54.6 


4.0 


41.4 


88.6 



Table 1 shows that multipole vectors associated to quadrupole, octupole and £ = 5 are 
mildly unlikely orthogonal to the dipole. Multipole vectors associated to £ = 4: are unlikely 
aligned to the dipole at the 99.1% of C.L.. 

In analogy to what performed in [5, 7], when 521 and S31 are combined as follows 

53"27 = (2 52i+3 53i)/5 (3.1) 

we find that the percentages to obtain a smaller value for the considered directions are 
0.8, 73.6, 2.6, 49.4 for DIP, NEP, NGP and NSGP respectively. This confirms that when 
quadrupole and octupole are considered together as in eq. (3.1) then the anomaly related 
to the dipole direction becomes more severe. Note that the estimator defined in cq. (3.1) 
differs from eq. (2.3) even if both are of course related in information content. Considering 
eq. (2.3), it is shown in [7] that the anomaly is as low as 0.3%. Therefore both estimators 
(i.e. S^2r <S'32i) provide the same level of anomaly^. 

Results for the second MC are reported in Figures 1, 2, 3, 4, 5, 6, 7 and 8 where we 
show the probability distribution function (pdf) for each Sn, for S'321 and for S^2i , and 
the value of the same estimator for the ILC map with various directions, represented by 
vertical lines in all the plots. We have chosen the following conventions: black line for the 
dipole, green line for north Ecliptic pole, blue line for the north Galactic pole, yellow line 
for the north Super Galactic pole. In Fig. 1 we recover the pdf of S'321, given in [7]. See 
appendix A.l for comments about its particular shape. In Fig. 2 we show the pdf for 521 

'^The Healpix parameter Nside is related to the number of pixels Npi^ through Npix = 12 N^ide- 

*Wc should remind that the ILC map used in [7] was the WMAP 3 yr ILC map and not the last one 

available now (analyzed in this paper). 



-4- 



5000 - 
4000 - 
3000 - 




0.0 0.2 0.4 0.6 0.8 1.0 

Figure 1: Copi et al.'s estimator [7], called 5'32i in this paper. The black vertical line represent 
the Dipole direction. The plot presents the counts (y-axis) versus the statistic (x-axis). 

4000 p ' ' ' I j ' ' I ' ' ' I ' ' ' I ' ' ' rz 




0.0 0,2 0.4 0.0 0.8 1.0 

Figure 2: 821- Vertical lines represent specific directions. We use the following conventions: black 
for the dipole, green for the north Ecliptic pole, blue for the north Galactic pole and yellow for the 
north Super Galactic pole. The plot presents the counts (y-axis) versus the statistic (x-axis). 



and in Fig. 3 for S31. Since the DIP-values stand in left tails of the pdfs of both estimators, 
it is possible to state that the multipole vectors associated to quadrupole and octupole are 
nearly orthogonal to the dipole. The probability to obtain a smaller value is 6.9% and 7.8% 
for and respectively. When these estimators are combined as in eq. (3.1) than we 
obtain the pdf shown in Fig. 8. In this joint case the probability to obtain a smaller value 
than the DIP value is 6.6%. In Fig. 4 it is shown the S'41 estimator. Since the DIP-value 
of this estimator is in the right tail of its pdf, the multipole vectors associated to £ = 4 
are unlikely aligned to the dipole. The probability to obtain a larger value is 0.9%. In 
Fig. 5 we show the ^51 estimator. The DIP-value for this estimator is in the left tail of the 
corresponding pdf. Therefore we conclude that the multipole vectors associated to ^ = 5 
are unlikely orthogonal to the dipole. The probability to obtain a smaller value is 0.3%. 
In Figs. 6 and 7 we show the pdfs for 5*61 and 871 respectively. In these cases, since the 
DIP-values fall in the middle of the two pdfs we do not find any anomalies. For sake of 
completeness, the probability to obtain a smaller value than the one that is observed is 
68.9% for the Sqi and 55.1% for the 5*71. In Table 2 the probabilities for each estimator 
and for each considered direction are summarized. The most unlikely percentages have 
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Figure 3: 5*31 . Same conventions as Fig. 2. Tlie plot presents the counts (y-axis) versus the 
statistic (x-axis). 




Figure 4: 5*41 . Same conventions as Fig. 2. The plot presents the counts (y-axis) versus the 
statistic (x-axis). 




Figure 5: 5*51. Same conventions as Fig. 2.. The plot presents the counts (y-axis) versus the 
statistic (x-axis). 



been reported in bold. 

In Table 3 we report the multipole vectors associated to the low multipoles of the ILC 
5yr map. For each multipole order ^ we give the i multipole vectors writing their colatitude 
and the longitude (in degrees). 
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0.0 0.2 0.4 0.6 0.8 1.0 

Figure 6: Sei. Same conventions as Fig. 2. Tlie plot presents the counts (y-axis) versus the 
statistic (x-axis). 
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Figure 7: 5*71 . Same conventions as Fig. 2. The plot presents the counts (y-axis) versus the 
statistic (x-axis). 




0.0 0.2 0.4 0.6 0.8 1.0 



Figure 8: 8^21 ■ Same conventions as Fig. 2. The plot presents the counts (y-axis) versus the 
statistic (x-axis). 

4. Conclusions 

Taking into account a new estimator defined in eq. (3.1), we have confirmed the unlikely 
alignment of quadrupole and octupole with the dipole in the WMAP Syr ILC map. More- 
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Table 2: Second MC: probabilities (in percentage) to obtain a smaller value for the considered 

directions and estimators 

Estimator DIP NEP NGP NSGP 



5'21 


6.9 


45.5 


9.6 


44.5 


S31 


7.8 


90.8 


12.4 


45.7 


anew 
'-'321 


6.6 


65.7 


10.5 


41.7 


S41 


99.1 


2.2 


55.6 


67.5 




0.3 


16.4 


69.0 


96.9 




68.9 


51.8 


37.7 


96.1 


Sri 


55.1 


11.3 


45.3 


79.4 



over we have shown that there are new directional anomahes. Multipole vectors associated 
to £ = 4 are unhkely ahgned with the dipolc at 99.1% C.L.. Note that this anomaly holds 
for both the MCs. Furthermore the multipole vectors associated to i = 5 are unlikely 
orthogonal to the dipole itself, with a probability at 99.7% C.L. when the MC is performed 
considering frozen the relative alignments of the low multipoles. This has been shown 
through MC simulations that properly take into account WMAP's anisotropic pixel noise 
level (see footnote 5). 

What causes these directional anomalies? It is difficult to answer this question. It is 
still unknown whether these anomalies come from fundamental physics or whether they are 
the residual of some not perfectly removed astrophysical foreground or systematic effect 
[13, 14, 15, 16, 17, 18, 19]. As an example of the latter kind, in references [20, 21] it 
is presented a stTidy about the impact of the dipole straylight contamination on the low 
amplitude of the quadrupole and on the low £ alignments for Planck ^ characteristics and 
capabilities. 

Of course these anomalies might come from fundamental physics, e.g. a non trivial 
topology of the universe [23] or magnetic fields [24]. See also [22] where a Quadratic 
Maximum Likelihood method has been adopted to study a dipolar modulation. 

As far as we know few attempts have been tried to study alignments anomalies in 
polarization [25] with WMAP data [26]. New data (i.e. Planck data) are awaited with 
great interest. 
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Table 3: Multipole vectors for the ILC 5yr map at low i 



Multipole order 


Colatitude ( j 


T J i'0\ 

Longitude [ ) 


1 = 2 


102.03 


305.16 


i = 2 


106.96 


184.22 


£ = 3 


77.98 


313.59 


^ = 3 


130.23 


272.81 


^ = 3 


97.18 


203.91 


£ = 4 


25.05 


203.82 


1 = 4: 


118.63 


154.69 


£ = A 


90.00 


68.91 


£ = 4 


128.68 


37.97 


i = 5 


57.20 


288.28 


£ = 5 


52.83 


39.38 


£ = 5 


125.69 


279.84 


£ = 5 


37.20 


228.60 


£ = 5 


92.39 


355.78 


£ = 6 


37.20 


34.20 


£ = 6 


54.32 


77.34 


£ = 6 


73.04 


288.28 


£ = 6 


81.61 


222.19 


£ = 6 


32.60 


239.32 


£ = Q 


105.71 


151.88 


£ — 7 


55 77 


■^45 Q4 


£ = 7 


69.26 


92.81 


£ = 7 


97.18 


203.91 


£ = 7 


48.19 


7.03 


£ = 7 


62.72 


174.78 


£ = 7 


157.94 


45.00 


£ = 7 


69.26 


286.88 
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A. Some technicalities about MC simulations 

In this Appendix we deal with some technicalities about the performed simulations. In 
particular we focus on the second aforementioned MC in which the relative alignments of 
the low multipoles are fixed and given by the ILC map. Subsections A.l and A. 2 refer to 
such a MC. 
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Figure 9: Bootstrapping on iS'321 estimator. The plot presents the counts (y-axis) versus the 
statistic (x-axis). 

A.l Shape of the estimators 

The shape of the pdfs of the estimators is due to the relative ahgnments of the multipole 
considered. In this Section we perform a sort of "bootstrap" to the ILC map in order not 
to change its angular power spectrum while obtaining a new map with different relative 
alignments among the multipole vectors. As an example, we recompute the pdf of 5321 for 
this new map. This is done just to show that the particular (and sometimes puzzling) shape 
of the pdf of the considered estimators strongly depend on the phases of map taken under 
analysis. In order to not to change the angular power spectrum, we arbitrarily reshuffle 
the aim for £ = 2, 3 of the ILC map as follows: 

n^LC ILC 

"2 1/2 ~^ "2 2/1 

and 

n^LC ILC 
"3 1/2/3 ~^ "2 2/3/1 

leaving a2^Q unchanged. In Fig. 9 we show how the shape of the pdf for the estimator 
5321 change. 

A. 2 Impact of the degradation procedure 

In this Section we briefly discuss our MC inputs. The ILC map is provided by the WMAP 
team at Ngide = 512 and has to be degraded to be handled in a MC. We chose to work 
at Nside = 16 and we have to keep under control the impact of this degradation^*^. In 
order to evaluate the effect of the degradation on our results, we have compared the pdfs 
of the considered estimators obtained performing a MC directly at high resolution (i.e. 
^side = 512) to the ones obtained by a (far less expensive) MC at low resolution (i.e. 
^side = 16). Since the high resolution MC is computationally heavy, we have reduced 

^"Since we are interested in quantities that depend on phases rather than norms of aim, the impact of the 
degradation upon angular power spectrum is not what we have to study. However, for sake of completeness, 
the effect of this degradation on angular power spectrum is of the order of few fiK^ that is well below the 
cosmic variance but above the angular power spectrum of the WMAP noise expected at those angular 
scales. 
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0.0 0.2 0.4 06 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 

Figure 10: Left Panel: Degradation effect on 5*21 . Right Panel: Noise effect on 5*21. Both the 
panels present the counts (y-axis) versus the statistic (x-axis). 



the number of random (angles) extractions to 10'^. As an example, in the left panel of 
Fig. 10 we show the difference between the histogram at high resolution and the histogram 
at low resolution for the estimator 5*21 . This plot of the difference has to be compared 
to the difference between the histogram at low resolution that include an (anisotropic) 
noise realization (at the V band WMAP level) and the histogram at low resolution (just 
degraded, without any noise realization on it). See right panel of Fig. 10 for the noise effect 
on the estimator S'21. The level of the fluctuations in the right panel of Fig. 10 is larger 
than in the left panel of the same Figure, showing that the effect of noise covers the effect 
of the degradation. In order to quantify degradation and noise effects, we compute the 
standard deviation of the two histograms of the differences, obtaining 0.353 for the noise 
impact, and 0.136 for the degradation impact. See Table 4 for a summary of the above 
defined standard deviation for all the considered estimators. Since the effect of V band 
WMAP (anisotropic) noise is larger than the effect of the degradation for all the considered 
cases, we conclude that the process of the degradation does not introduce any significant 
spurious effect on the estimators we used to study anomalous directions. 



Table 4: Impact of noise vs degradation 



Estimator 


^ Noise 


^Degrad 


5'21 


0.353 


0.136 


S31 


0.444 


0.066 


S41 


0.403 


0.108 


S51 


0.242 


0.107 




0.573 


0.191 




0.512 


0.132 



References 



[1] E. Komatsu et al. [WMAP CoUaboration] , arXiv:0803.0547 [astro-ph]. J. Dunkley et al. 
[WMAP Collaboration], arXiv:0803.0586 [astro-ph]. 



[2] C. J. Copi, D. Hutcrcr and G. D. Starkman, Phys. Rev. D 70, 043515 (2004) 

[arXiv:astro-ph/0310511]. 

J. R. Weeks, arXiv:astro-ph/0412231. 

M. Tegmark, A. de Oliveira-Costa and A. Hamilton, Phys. Rev. D 68, 123523 (2003) 
[arXiv:astro-ph/0302496]. D. J. Schwarz, G. D. Starkman, D. Huterer and C. J. Copi, Phys. 
Rev. Lett. 93, 221301 (2004) [arXiv:astro-ph/0403353]. K. Land and J. Magueijo, Phys. Rev. 
Lett. 95, 071301 (2005) [arXiv:astro-ph/0502237]. C. Vale, arXiv:astro-ph/0509039. 

C. J. Copi, D. Huterer, D. J. Schwarz and G. D. Starkman, Mon. Not. Roy. Astron. Soc. 
367, 79 (2006) [arXiv:astro-ph/0508047]. 

A. Gruppuso and C. Burigana, JCAP 0908, 004 (2009) [arXiv:0907.1949 [astro-ph.CO]]. 

C. Copi, D. Hutcrcr, D. Schwarz and G. Starkman, Phys. Rev. D 75, 023507 (2007) 

[arXiv:astro-ph/0605135]. 

L. R. Abramo, A. Bernui, I. S. Ferreira, T. Villela and C. A. Wuensche, Phys. Rev. D 74, 
063506 (2006) [arXiv:astro-ph/0604346]. 

G. Katz and J. Weeks, Phys. Rev. D 70, 063527 (2004) [arXiv:astro-ph/0405631]. 

K.M. Gorski, E. Hivon, A.J. Banday, B.D. Wandelt, F.K. Hansen, M. Reinecke, and M. 
Bartelmann, HEALPix: A Framework for High-resolution Discretization and Fast Analysis of 
Data Distributed on the Sphere, Ap.J., 622, 759-771, 2005. 

K. Land and J. Magueijo, Mon. Not. Roy. Astron. Soc. 362, 838 (2005) 
[arXiv:astro-ph/0502574] . 

M. R. Dennis, J. Phys. A 37 (2004) 9487-9500, J.Phys. A38 (2005) 1653-1658 

L. R. Abramo, L. S. Jr. and C. A. Wuensche, Phys. Rev. D 74, 083515 (2006) 

[arXiv:astro-ph/0605269]. 

P. D. Nasclsky and O. V. Verkhodanov, Int. J. Mod. Phys. D 17, 179 (2008) 
[arXiv:astro-ph/0609409] . 

K. T. Inouc and J. Silk, Astrophys. J. 648, 23 (2006) [arXiv:astro-ph/0602478]. 
K. T. Inouc and J. Silk, Astrophys. J. 664, 650 (2007) [arXiv:astro-ph/0612347]. 
A. Cooray and N. Seto, JCAP 0512, 004 (2005) [arXiv:astro-ph/0510137]. 

H. V. Peiris and T. L. Smith, arXiv: 1002.0836 [astro-ph.CO]. 
C. L. Francis and J. A. Peacock, arXiv:0909.2495 [astro-ph.CO]. 

C. Burigana, A. Gruppuso and F. FinelU, Mon. Not. Roy. Astron. Soc. 371, 1570 (2006) 

[arXiv:astro-ph/0607506] . 

A. Gruppuso, C. Burigana and F. Finelh, Mon. Not. Roy. Astron. Soc. 376 (2007) 907 

[arXiv:astro-ph/0701295] . 

D. Hanson and A. Lewis, Phys. Rev. D 80, 063004 (2009) [arXiv:0908.0963 [astro-ph.CO]]. 
P. Bielewicz and A. Riazuelo, arXiv:0804.2437 [astro-ph]. 

A. Bernui and W. S. Hipolito-Ricaldi, Mon. Not. Roy. Astron. Soc. 389, 1453 (2008) 
[arXiv:0807.1076 [astro-ph]]. 



- 12 - 



[25] C. Dvorkin, H. V. Peiris and W. Hu, Phys. Rev. D 77 (2008) 063008 [arXiv:071 1.2321 

[astro-ph]]. 

[26] M. Frommert and T. A. Ensslin, arXiv:0908.0453 [astro-ph. CO]. 



-13- 



